March 1998 



UND-HEP-97-US03 
FTUAM 98/3 



hep-ph/980345Tl 



Imminent Phenomenology of a Minimal Gauge-Mediated Model 



00 

[ Emidio Gabrielli* 

G\ _ 

Departamento de Fisica Teorica 
Universidad Autonoma de Madrid 
Cantoblanco, 28049 Madrid 
Spain 

Uri SaricU 

> . 

Department of Physics 

University of Notre Dame 

^ ■ Notre Dame, IN 46556 

O : USA 
00 

^ ' Abstract 

a: 

Q_i| We calculate the inclusive branching ratio for B — > X s j, the inclusive branching 

ratios and asymmetries for B — > X s £ + £~, and the anomalous magnetic moment — 2 
of the muon, within a minimal gauge-mediated SUSY-breaking model which naturally 
generates a large ratio tan (3 of Higgs field vacuum expectation values. These predictions 
are highly correlated with each other, depending on only two fundamental parameters: 
the superpartner mass scale and the logarithm of a common messenger mass. The 
predictions for B — > X s j decay and — 2 are in somewhat better agreement with current 
experiments than the standard model, but a much sharper comparison will soon be 
possible using new measurements now in progress or under analysis. Moreover we predict 
large deviations in B — > X s e + e~ and B — > X s /j, + /i~ asymmetries, and somewhat smaller 
ones in B — > X s e + e~ and B — > X s t + t~ branching ratios, which will be detectable in 
hadronic colliders. 
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1 Introduction 



Arguably the least-understood aspects of supersymmetry (SUSY) as a physically-acceptable 
theory are the mechanisms which break SUSY and which communicate that breaking to the 
observable sector. These mechanisms must generate soft-breaking mass terms which are small 
enough to protect the Higgs mass, large enough to have evaded all direct searches so far, 
and flavor-symmetric enough to respect all current indirect bounds. In particular, flavour- 
changing neutral current (FCNC) processes impose strong constraints on any non-universal 
soft-breaking masses, due to gluino- and photino- mediated SUSY contributions Moreover 
CP-violating phases in the soft-breaking sector must be strongly suppressed to satisfy the 
constraints on electric dipole moments 

In supergravity-mediated scenarios with flat Kahler metric H the SUSY breaking is com- 
municated from the hidden sector to the observable one by gravitational interactions, leading 
to universality in the soft-breaking sector near the Planck scale. But even such a strong as- 
sumption about physics at extremely high scales may be not sufficient to prevent large flavor 
violations at observable scales. Indeed in grand-unified supersymmetric theories , the large 
top Yukawa couplings and radiative effects above the unification scale can generate unaccept- 
ably large FCNC and CP-violating processes. Moreover it is doubtful whether supergravity 
theories derived from superstrings would have flat Kahler metrics and flavour-independent 
supersymmetry-breaking terms ||. 

Recently there has been a revival of interest in a class of models which can solve both the 
SUSY flavour and CP problems: theories of gauge-mediated SUSY breaking 0-f|. In these 
theories the SUSY-breaking sector is coupled to the observable (minimal supersymmetric stan- 
dard model, or MSSM) sector via gauge interactions. Many models employ the standard-model 
(SM) gauge interactions to couple the MSSM to a messenger sector, which then couples more 
directly to the SUSY-breaking sector (much as in the gravitationally-mediated picture the 
MSSM and the SUSY-breaking hidden sector communicated only via Planck-scale physics). 
The messenger sector may be far below the Planck scale, thus also lowering the SUSY-breaking 
scale, and the gauge interactions which generate observable SUSY-breaking masses are natu- 
rally flavor-blind, thus preserving the approximate flavor symmetries of the SM. In the class 
of models studied in this work, the ultimate source of SUSY breaking is parametrized by a 
superfield S whose scalar and F components both acquire vacuum expectation values (VEVs) 
through an unspecified mechanism. Since this superfield is a SM gauge singlet with no direct 
couplings to the observable sector, the mechanism of SUSY breaking is largely isolated from, 
and therefore unconstrained by, the known properties of the SM. S does couple to a set of 
messenger fields, which do carry SM gauge quantum numbers. When their superpartners are 
split by the couplings to S, they radiatively split the MSSM gauginos from the gauge bosons 
and the sfermions from the SM fermions, thus conveying SUSY breaking to the particles we 
observe. 

The pattern of SUSY breaking in the MSSM produced by gauge mediation is thus exactly 
what is needed to solve the SUSY flavor and CP problems. Moreover the most minimal 
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of these models are highly predictive, since most of the MSSM parameters beyond the SM 
depend on only the scale of the messenger masses and of the SUSY-breaking VEV (and the 
number of messenger fields). Unfortunately neither the superpotential /i term, which is the 
coupling between the up- and down-type Higgs superfields, nor the soft-breaking term B/i, 
which couples the scalar Higgs fields, is generated by the usual gauge mediation mechanism: 
both couplings violate the Peccei-Quinn (PQ) symmetry of relative rotations between the 
Higgs doublets, which gauge interactions preserve. But for phenomenological reasons /i must 
be generated somehow, and indeed a few dynamical mechanisms have been proposed to break 
the PQ symmetry at the proper scale ||,||. Possible solutions of this \i problem are an active 
area of research in this field. 

We will regard the \x term as a free parameter, to be fixed phenomenologically by requiring 
that the electroweak vacuum be broken at the correct scale by radiative corrections. We will 



then assume, as in ref. ||10|| , that when the /x term is generated the corresponding B/j term is 
not (and also the soft-breaking Higgs masses are unaltered). Of course B\i will be induced 
radiatively, but since we assume that its boundary value at the messenger scale vanishes, we 
can predict its value at lower, observable scales ||10|| . The vanishing of B/j, at the messenger 
scale offers several attractions: 



1. The two physical SUSY CP violating phases axg(BM*/ 2 ) and arg(AM*/ 2 ), where M1/2 is 
a gaugino mass and A is the scalar analog of the Yukawa couplings in the superpotential, 
are zero @, solving completely the SUSY CP problem; 

2. all soft-breaking terms and relative signs can be predicted essentially in terms of only 
two parameters, a common messenger mass and the SUSY-breaking scale; 

3. Since Bfi at observable scales turns out to be small, tan j3, the ratio between the up- and 
down-type Higgs VEVs, is naturally predicted to be large, thus generating the observed 
large hierarchy between top and bottom quark masses without appealing to hierarchical 
Yukawa couplings. 



Since a naturally large tan (3 is a signature of this model, we analyse a class of low energy 
processes which are particularly sensitive to large tan/3. It is well known that in the MSSM 
the processes mediated by magnetic dipole transitions can be enhanced by tan/?. The cor- 
responding SM amplitudes are chiral suppressed, while some of the SUSY contributions can 
receive a tan/3 enhancement, making the new contributions competitive with the standard 
ones. Within this class of processes, we consider in the framework of the MGM model the 
radiative b —>■ 57 and semileptonic b — ► sl + l~ decays (with I = e, /i, r) and the contribution to 
the anomalous magnetic moment of the muon — For the b —> s£ + £~ processes we analyse 
both the decay rate and the energy or forward-backward (FB) asymmetries. 

In b —>■ 57 and — 2 the corresponding amplitudes are proportional to the magnetic 
dipole transitions and so are directly enhanced by tan/3. The b — ► sl + l~ amplitude receives 
contributions from both the local four-quark operators Qg and Qio, which are not enhanced by 
tan /3 and in fact will be largely unaffected by the MGM, and from the same magnetic-dipole 
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operator Q 7 responsible for the b — > sj decay. (Qj contributes to the b — > s£ + £~ amplitude 
through the exchange diagram between Q 7 and the electromagnetic lepton current.) Therefore 
the sensitivity of b — ► s£ + £~ to tan f3 is correlated with that of b —>■ 37 through their common 
dependence on a single Wilson coefficient. 



A number of studies for b -> 57 |P | -[g7 | , b -> s£ + £~ 111 -PI 



P and ^ - 2 g|-| 

have already been carried out for MSSM models similar to the one we consider. But our 
analysis for the MGM model is unique in its predictivity: there are essentially only two 
parameters, and they determine a tight correlation between these processes. (Moreover in our 
analysis we include the complete next-to-leading order accuracy in the strong coupling for 

a 



the SM contribution to b 



S7 



and b -> sl + l~ |27 



Present data mildly favor our 



predictions over those of the SM. Measuring the correlated deviations in these processes from 
SM predictions will either favor this model over the standard one and measure its parameters, 
or strongly constrain the model by forcing it to agree closely with the SM, or disfavor both 
this model and the standard one. 

The paper is organized as follows : in the next section we define the model and give 
the analytical results for the mass spectrum of squarks, sleptons, Higgs and gauginos at the 
physical scale. We state these results for the mildly more general model than the minimal 
one, allowing for N > 1 messenger families. All these quantities are expressed as a function of 
the SUSY-breaking scale A (or equivalently the wino mass) and the common messenger mass. 
Also in this section we give the numerical predictions for the magnitude of the \x parameter 
(for various values of N, messenger scale and top mass) and a detailed discussion about the 
predictions of its sign. In Sees. 3, 4 and 5 we analyse the contribution of the MGM model to 
the b — > sj, b — > s£ + £~ decay rates and asymmetries, and — 2, respectively. We give model- 
independent parametrizations for the total rates and asymmetries and analytical expressions 
for the MGM contributions to the corresponding amplitudes in a suitable approximation. In 
Sec. 6 we present and analyse the numerical predictions for the b — > S7 and b — > s£ + £~ decay 
rates and asymmetries versus the — 2 ones, as a function of the SUSY breaking scale and 
the common messenger mass in both N — 1 and N = 2 scenarios. The final section contains 
our conclusions and outlook. 



2 Minimal gauge mediated model 

In GMSB models supersymmetry is broken at some scale above the electroweak scale and then 
the breaking is communicated to the MSSM particles by the usual SM gauge interactions |J- 
||. In the minimal realization of this idea, a pair of SU(2) doublet chiral superfields $2; $2 
and a pair of SU(3) triplet chiral superfields $3, $3, called messenger fields, couple to a singlet 
chiral superfield S with Yukawa couplings A2,3- The messenger fields get their supersymmetric 
mass through the VEV (S) of the scalar component of the S superfield. But the F component 
of S also has a nonvanishing VEV, (F), which parametrizes the breaking of supersymmetry. 
The breaking is manifest in the consequent mass splitting between the fermionic and scalar 
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component of the messenger fields the fermions of $; acquire masses Mm, = Xi(S) while 
their scalar partners get squared masses mf = \\i(S)\ 2 ± |Aj(F)|. 

In this minimal gauge-mediation (MGM) model, the messenger fields couple at tree level 
only with the gauge sector of the MSSM and not with its matter fields. Thus SUSY breaking 
is transmitted to the MSSM only by the ordinary gauge interactions. At one loop, gaugino 
masses are generated by finite self-energy diagrams, with the fermionic and scalar components 
of the messenger fields running in the loop. In the following we will allow the possibility of 
N > 1 messenger families (pairs of doublets and triplets), but only a single source S of SUSY 
breaking. We will also assume that the messengers are approximately degenerate; our results 
are not very sensitive to this splitting as long as the messengers are not very split and not 
very close to their lower bound (see below). We denote quantities evaluated at the messenger 
scale by an overbar. Then the gaugino mass spectrum M$ is given by 

Mi = N-±Kg{ Xi ) = M l9 (x t ) (1) 

47T 

where A = (Fs) / (S) and Xi = K/Mm v The function g(x) contains the result of the one-loop 
integration. Its argument must satisfy Xi < 1 to ensure that the messenger scalars have a 
positive mass-squared. We will not be concerned with Mm, just above but very close to the 
minimal value A, so Xi will never approach unity, in which case g{xi) ~ 1 and ~ Mi will 
be very good approximations. 

At two-loop order in the messenger fields, or more plainly at one-loop order in the gaugino 
masses, squark, slepton and Higgs SUSY-breaking squared masses are generated by self-energy 
diagrams with a gaugino mass insertion. Thus the sfermion masses (rather than their squares) 
are the same order as the gaugino masses. Once again we are concerned with X{ well below 
unity, in which case the results of the sfermion mass calculation can be expressed as [ ID ]: 



ml ~ i (2C3MI + 2C 2 M* + ~ Y 2 M 2 ) (2) 

where the coefficients Cj depend on the gauge quantum numbers of the scalar multiplets and 
Y is the weak hypercharge. Explicitly, these equations read: 

r?? = 1 — (-Ml + -Mf + —Adf) (3) 
Q N \3 2 30 / V ' 

™l = ll(~Mi + ^M?) (4) 



N \3 15 

— -Ml H . 

N V3 3 15 

— -Ml + —. 
N \2 2 10 



™I = 1 -r(oMl + -Ml) (5) 
K = ^(S^ + ^Mi) (6) 



^2 , 1 {6nr2 



= 1-t-Mf) (7) 

m 2 H =m 2 H = -(-MI + —M 2 ) (8) 
Hu h d N \2 10 J K J 
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for the squared masses of the squark doublet, the up- and down-type squark singlets, the 
slepton doublet, the charged slepton singlet, and the up- and down-type Higgs doublets. 
We've also explicitly inserted the unit matrix 1 in flavor space: the gauge interactions which 
generate these masses are flavor-blind. 

Another consequence of the pure gauge mediation in the MGM is that the SUSY-breaking 
A terms which couple the squarks (or sleptons) to the Higgs doublets are not induced at the 
same order as the gaugino or sfermion masses, since they break the chiral (flavor) symmetries 
while the gauge interactions do not. A terms are induced at higher order via gaugino masses, 
so both A and M arise through the same source A, and hence the relative CP-violating phase 
between them axg{A*M^) vanishes naturally. 

The \x parameter, which couples the up- and down-type Higgs doublet superfields in the 
superpotential, is not generated within the GMSB paradigm, since it breaks the Peccei-Quinn 
symmetry while the gauge interactions do not. To have a (phenomenologically- mandated) 
nonvanishing /i term, the MGM model must be extended. We will assume that such an 
extension does not have any impact on the mass parameters beyond the // term (but see Ref. || 
for counterexamples). Thus effectively [i is an arbitrary additional parameter of the model, 
as was assumed in most of the previous literature. But in fact its value can be predicted by 
the phenomenological requirement that the electroweak vacuum should spontaneously break 
to yield the proper Higgs VEV. 

The SUSY-breaking analogue of the /z term, the Bfi term coupling the Higgs doublets in the 
scalar potential, is also not generated by the pure MGM. We will assume, as we did above, that 
whatever mechanism generates \x does not generate B/j,, and this is the major assumption of 
this work. While this is by no means a generic aspect of GMSB models (see Ref. |§ for critical 
analysis of this issue), it is theoretically conceivable since Bfi would break an R symmetry 
which /i would not (see Ref. |7| for an example and Ref. |10| for further discussion) . It is also 
phenomenologically appealing, since (as for the A terms) B is generated radiatively through 
gaugino masses, and hence the CP-violating relative phase axg(B* M3) once again vanishes 
naturally. Thus the SUSY CP problem is completely absent in this model. To reiterate: in 
the following we will always assume the boundary condition B ~ 0. 

Notice that the only source of phases (apart from topological vacuum angles) is the SM 
CP- violating phase in the quark Yukawa couplings. In particular, no physical relative signs 
are arbitrary parameters; thus the relative sign between fi, the Higgs VEVs and the gaugino 
masses (which is often just called the sign of /i) is a prediction of this model, which we calculate 
below. 

Since the messenger scale is considerably higher than the electroweak scale, we use the 
above MSSM parameters as boundary conditions to the renormalization group (RG) evolution 
equations, and calculate the RG-improved masses closer to the electroweak scale. But as long 
as the messenger scale is only a few orders of magnitude above the electroweak scale, a one- 
step solution to the RG equations will suffice for all parameters except B, for which a delicate 
cancellation occurs. The RG evolution starts at the messenger scale; it effectively stops at a 



typical squark mass scale m. Both of these statements were discussed and justified in Ref. |10[; 
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in particular, the doublet and triplet messenger masses are expected to be quite similar in 
a unified model, so there is a reasonably well-defined single messenger scale, and most of 
the results do not depend significantly on the small splitting between Mm 2 an d Mm 3 . We 
will usually only keep the leading powers of L/(8n 2 ) = ln(M M /m)/(87r 2 ), where M M is the 
messenger mass scale. Since tan (3 = (Hy)/{H D ) will also turn out large, we will also usually 
keep only leading powers of 1/ tan/3. 

The gaugino masses, and thence the physical chargino and neutralino mass matrices, are 
given at the squark scale respectively by: 



Mi = 



M x± 



M x o 



^)Mi 



M 2 V2M W 
fi 
/ M 1 







M 2 




\sin 9wmz — cos Owmz 







-fj, 



sm6 w m z ^ 
— cos 6 w m z 







/ 



(9) 
(10) 

(11) 



The up- and down-type squark mass matrices M~ and Mj respectively at the squark scale 



are 



where 



Ml = [ m Q 



mi 



§ sin 2 



6w) 



-m u A v 

m| + m 2 , + 1 (| - \ sin 2 W 
lirud tan (3 



to| 



mi 



TO- 



TO 



TO, 



-m u Ay 

2 i 2 „;„2 



l^sin 9 w m z 



/iTO^tan p 



m\ 



l„-sin Owf^z 



(12) 
(13) 



TOa 



Q 
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1^2 1 o^r^ 



= to^ + \ 1 ( -^M; + 3^M; + -^M; ) - ^ ( -M 3 ^ + 3M^ + — Mf 



X 2 
N 



16 - 9 - n 7 « o X 

-Ml + 3M 2 + -M 2 ^ 



15 
L 

8^ 



A 2 /16 



13 



15 



(14) 



ml = ml + \1 (fgjMl + ^ 2 M?) - ^ (^M 3 2 + 3M 2 + ^M 2 ) } A (15 ) 



n + U ( + ±^M?) - ^ (^M 3 2 + 3M 2 2 + ^ 



8tt 2 ' 



(16) 



to m (A;y) and mj (Ad) are the up- and down-type quark mass (Yukawa coupling) matrices in 
flavour space, and Ay is given by 



(16 13 \ L 

—gl~M 3 + 3glM 2 + —gfM 



8tt 2 



(17) 
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where the boundary condition Ajj ~ has been used. 

For the slepton mass matrices we find, with the same approximations as for the squarks, 



where 



Ml 



Mi 



Tflr 



m-, 



+ m 2 + 1 f | — sin 2 6 W J 
firni tan (5 

— l-m% 



111; 



mf + li 



1 ( 3g 2 2 M 2 2 + 1<£m[ 



12 n ^^2\ 2A 2 



-9iMl - 



fimi tan (3 
mf + 1 sin 2 Ow'm'z 



N 



3M2+ ^ 2 )}^ 



(18) 
(19) 

(20) 
(21) 



and rrii {^i) are the diagonal lepton mass (Yukawa) matrix in lepton flavour space. 

The remaining parameters to be determined are the Higgs couplings fi and Bfi. These 
enter the scalar potential and determine the VEVs of the two Higgs doublets vu and in 
terms of the soft-breaking masses (discussed above) and the quartic Higgs couplings (fixed 
supersymmetrically by the gauge couplings). Since we are assuming B = 0, the only free 
parameter is /i, and it is set by requiring agreement with the experimental observable v 2 = 
vfj+vjj ~ (174 GeV) 2 , or equivalently with m| = (gf+g^v 2 /2. Since the coupling B\i between 
the up- and down- type Higgs scalars will be very small, v ~ vu ^> v^, so the electroweak scale 

requirement is essentially a condition on the up-type Higgs mass-squared: m 2 Hu = — m 2 z 

to leading order. Substituting in the up- type Higgs mass as a function of /z 2 , and including 
higher-order corrections, leads to: 



/i 



A^A 2 



m 



H 



1 

-m 
2 



1(1 + ^ 



(22) 



1) accounts for the loop correc- 

— + -M 2 ) + 

N\2 2 10 V 



The details of this minimization are given in Ref. [T^j: 5h 
tions to the Higgs quartic couplings induced by the heavy top, m 2 H 

?>g\M\ (^—^j is the common Higgs mass, and A 2 A 2 is the Yukawa coupling correction which 

lowers the up-type Higgs mass and breaks the electroweak symmetry. Neglecting terms smaller 
than ~ 1%, and recalling that X t is evaluated at the squark mass scale as in Ref. ]I0 |, we find: 



A 2 - 



3_La 

N8n 2 



16 



M\ + 3M 2 



+ 



L. 



in' 



NfgWl - (™9l + 8gj ) M 3 2 - 8^ 2 M 2 



(23) 



We will set the logarithm to be La = ln(MM 3 /rriQ) +3/2 where the constant term, calculated 
in Ref. ||13|| , accounts for the 1-loop thresholds at the messenger scale. The absence of the 
complete 1-loop threshold expression resulted in the major uncertainty in the previous analysis 
of this model |K| . 
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71 /T /A 

%/A 


m 4 = loo (jrev 


rrit = 175 Uev 


rrit = 185 GeV 


1 


23 GeV + 1.08M 2 


23 GeV + 1.22M 2 


23 GeV + 1.36M 2 


2 


22 GeV + 1.46M 2 


23 GeV + 1.64M 2 


23 GeV + 1.82M 2 


10 


22 GeV + 1.62M 2 


22 GeV + 1.81M 2 


22 GeV + 2.01M 2 


100 


20 GeV + 1.70M 2 


20 GeV + 1.90M 2 


20 GeV + 2.11M 2 


1000 


18 GeV + 1.74M 2 


18 GeV + 1.95M 2 


17GeV + 2.16M 2 


10000 


15 GeV + 1.76M 2 


15 GeV + 1.97M 2 


14 GeV + 2.19M 2 



Table 1: The MGM predictions for /i as a function of the wino mass M 2 in the linear approx- 
imation and for some values of messenger mass Mm 2 and top mass m t , with Mm 3 = 1.3Mm 2 
and N = 1. 



M M JA 


m t = 165 GeV 


m t = 175 GeV 


m t = 185 GeV 


1 


22 GeV + 1.26M 2 


22 GeV + 1.43M 2 


22 GeV + 1.60M 2 


2 


21 GeV + 1.43M 2 


21 GeV + 1.61M 2 


21 GeV + I.8OM2 


10 


21 GeV + 1.57M 2 


21 GeV + 1.76M 2 


20 GeV+ 1.96M 2 


100 


19 GeV + 1.66M 2 


19 GeV + 1.86M 2 


19 GeV + 2.07M 2 


1000 


17 GeV + 1.71M 2 


17 GeV + 1.91M 2 


16GeV + 2.12M 2 


10000 


15 GeV + 1.73M 2 


15 GeV + 1.94M 2 


14 GeV + 2.15M 2 



Table 2: The MGM predictions for /i as a function of the wino mass M 2 in the linear approx- 
imation and for some values of messenger mass Mm 3 and top mass mt, with M\i_, = 1.3Mm 3 
and N — 1. 



Numerically, we find that /i increases approximately linearly with the superpartner scale, 
which we characterize by the wino mass M 2 , even for light M 2 . Specifically, we find that 
/i is given by the linear expressions of Tables [l|-44| for various values of N, the messenger 
scale and the top mass. These approximations yield fi within ~ ±10 GeV over the range 
100 GeV < M 2 < 400 GeV. Several features are evident in these results: /i increases with 
the messenger scale and with the top mass, since it must cancel a larger radiative correction 
; it decreases with TV relative to M 2 , since more messengers lower the sfermions relative to 
the gauginos [see Eqs. OE0)]; and it's rather insensitive to small messenger splittings induced 
presumably by RG evolution from high scales. 

Finally, tan/? = vjj/vr, is predicted by minimizing the full scalar potential for both Higgs 
doublets: 

— = ~^ (24) 
tan/? m 2 eSA 

A striking aspect of the MGM is that, assuming B = (or practically zero), the low-energy 
value of B is also very small, and hence tan f3 is naturally very large 0, [Kj. The reason is not 
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M M Jk 


m t = 165 GeV 


m t = 175 GeV 


m t = 185 GeV 


1 


8GeV + 0.81M 2 


9 GeV + 0.91M 2 


10 GeV + I.OIM2 


2 


8GeV+ 1.10M 2 


9 GeV + 1.23M 2 


10 GeV + 1.36M 2 


10 


9GeV+ 1.23M 2 


10 GeV + 1.37M 2 


11 GeV + 1.51M 2 


100 


9GeV+ I.3IM2 


10 GeV + 1.46M 2 


10 GeV + I.6OM2 


1000 


8GeV+ 1.35M 2 


9 GeV + 1.50M 2 


9GeV+ 1.65M 2 


10000 


7GeV+ 1.37M 2 


8GeV+ 1.53M 2 


8GeV+ I.68M2 



Table 3: The ji MGM predictions as in Table 1 with N = 2 



AW A 


mt = 165 GeV 


m t = 175 GeV 


m t = 185 GeV 


1 


7GcV + 0.96M 2 


9 GeV + I.O8M2 


10 GeV + I.2OM2 


2 


6 GeV + I.O8M2 


7GeV+ I.2IM2 


8GeV + 1.34M 2 


10 


7GeV+ 1.19M 2 


9 GeV + 1.33M 2 


9GeV+ 1.47M 2 


100 


8GeV+ 1.27M 2 


9 GeV + 1.42M 2 


9GeV+ 1.56M 2 


1000 


7GeV+ 1.32M 2 


8 GeV + \A1M 2 


8GeV+ 1.62M 2 


10000 


6 GeV + 1.34M 2 


7GeV+ 1.50M 2 


7GeV+ 1.65M 2 



Table 4: The \i MGM predictions as in Table 2 for iV = 2 
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simply a small amount of RG evolution between the messenger and squark mass scales, but 
rather jlO| a cancellation between the direct gaugino contribution to the RG evolution and the 



direct v4-term (or indirect gluino) contribution to that evolution — a fortuitous cancellation 
that happens when Mm is not too far above its minimal value A. Thus if the messengers are 
not too heavy, the B parameter evaluated at the squark mass scale is very small, as it is at 
the messenger scale, although for intermediate scales B is considerably larger. So the coupling 
between the Higgs scalars is predicted to be very small at the relevant scale of electroweak 
breaking, which leads to a large hierarchy between their VEVs: tan (3 — Vu/vd ^> 1- Because 
of the cancellation of the leading-order terms, the _B(mg) must be calculated quite carefully at 
the next-to-leading order. The analytic results of this calculation, and the expression for the 
effective pseudoscalar mass-squared parameter A to be used in minimizing the potential, 
are given in Ref. [|IU|| . Then tan/3 is predicted by solving Eq. ( |24| ) numerically, taking into 
account the implicit dependence of the right-hand side (and in particular of B) on tan /3. 

In addition to new mass scale parameters, the MSSM is characterized by two new phase 
parameters beyond the SM. The MGM predicts these signs in any fixed convention. To see 
this, note that they may be put entirely into the A and B parameters by appropriate field 
redefinitions ||12|| . But in minimal gauge mediation the A and B parameters of the effective 
MSSM at the messenger scale vanish, and hence the phase parameters are predicted: they 
vanish. Then their values at lower scales follow from the calculable RG evolution of A and 
B. In a convention where the gaugino masses and top and bottom Yukawa couplings are real, 
no complex phases can then be generated (at least in low orders) in A and B, only a relative 
sign, which we calculate below. 

This sign is particularly crucial for the observables of interest to us, g M — 2 and b — > 57, 
which are enhanced by large tan {3. As shown in the first reference of |T3j , all such observables 
proportional to tan (3 are also proportional to the degree of breaking of the Peccei-Quinn (PQ) 
symmetry of the Higgs sector, characterized by /z, and of an R symmetry, characterized by the 
gaugino mass M1/2. The amplitudes of these observables will be determined by the sign of the 
product tan/3/iM!/ 2 , which is a prediction of the MGM in any fixed Lagrangian convention, 
as we calculate below. The final predictions of the observables will of course be independent 
of any convention. 

To state the sign, we must establish some convention. In 2-component Weyl notation, we 
use a Lagrangian schematically of the form: 



£ gaug e ~ W [Gt L t L - Gt R t R J (25) 

^Yukawa / dPOW, W ~ fiH°H° D + \ t H°t L t R (26) 

Aoft M 1/2 GG - \ t A t t L t R H° - BtiH^H° D (27) 

C F ~ - \dW/d<P\ 2 ~ -X t fiH° D t L t R (28) 

Here Hy D are the neutral components of the Higgs doublets, G is the gaugino field, and 
the proper group structure and indices are implied: our purpose is only to establish the sign 
convention. With this Lagrangian, we calculate (and find agreement with Moroi j|7|]) that, 
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for example, the dominant MSSM contribution to — 2 is Aa ~ +Mi/2/xtan/3. We can also 
easily extract and minimize the scalar Higgs potential to find 1/ tan/5 ~ —Bfi. Next, a 1-loop 
top-gluino diagram using this Lagrangian generates A t ~ +M\/2- Then B is generated by 
both a 1-loop stop diagram with an A insertion, Ba ~ —A t ~ — M1/2, and a 1-loop higgsino- 
wino diagram, 5^ ~ +My2- The coincidental near-cancellation of these two is the reason 
tan/? is large in the MGM, but we find that essentially always \Bjs\ > \Bg\, so B ~ — M\/2- 
Therefore with this sign convention M^/^tan /3 > 0. Hence the convention- independent, 
physical sign of Aa is positive, and similarly the deviations of the various b — > S7 amplitudes 
are unambiguously predicted. 

To be more specific, for the remainder of this paper we will further specify that fi and 
M1/2 are positive, and therefore so is tan/3. To translate our intermediate results to different 
Lagrangian conventions, the interested reader need only use the Lagrangian of Eqs. fl25| ) - 
( ^81) to determine the desired signs of [i and Mi/ 2 , and then establish the corresponding sign of 
tan/? from M!/ 2 /xtan/3 > 0. Our final results are physical, and therefore will be unchanged. 
In much of the existing literature, the signs of M ly / 2 and tan/3 are fixed by convention, in 
which case the MGM makes a prediction of the sign of /z. 

We conclude this section by observing that the low energy mass spectrum in the MSSM, 
as well as tan/3 and the sign of /1, can be predicted within the MGM in terms of only the 
SUSY-breaking scale A (or equivalently any particular superpartner mass, say the wino M2) 
and the logarithm of the messenger mass scale Mm- This extraordinary predictivity forms the 
basis for the remainder of this work: the formulation of low-energy experimental tests of the 
MGM. 



3 b -> S7 

As discussed in the previous section, an interesting feature of the MGM model is that tan j3 
is predicted large. Since we are interested in constraining this model from low energy physics, 
we will restrict our investigation to some processes which are particularly sensitive to large 
tan/3 values. It is well known that the low energy processes mediated by magnetic dipole 
transitions can be significantly enhanced in the MSSM for tan/3 3> 1, because the chiral 
suppression of superpartner-mediated diagrams is removed to a large extent by the large 
tan/? (i.e. the Yukawa couplings are no longer small), which allows for a competition with 
the chirally-suppressed SM amplitude. Among this class of processes we consider here the 
rare B — > X s 7 decay. In particular in this section we will give the analytical results for the 
dominant contributions to the B — > X s ~f decay in the MGM model, and fold them in with the 
recent Next to Leading Order (NLO) calculation for the SM. 

Let us start with the experimental results. The most recent published result by the CLEO 
collaboration for the total inclusive branching ratio BK(B — > X s '-f) is p4f 

10 4 BR expt (£ _^ x s7 ) = 2.32 ± 0.67 (29) 
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combining 1 a uncertainties, and 

1.0 < 10 4 BR expt (S -> X s7 ) < 4.2 (30) 

at the 95% confidence level. 

Recently the NLO corrections to the quark b — > S7 decay rate have been completed in 
Ref. [^1[], including the calculation of the three- loop anomalous dimension matrix of the effec- 
tive theory. This is a necessary ingredient for a consistent resummation, at the NLO accuracy, 
of the large leading logarithms, log M w /m h . 

The inclusion of the NLO corrections in the h — > sj decay rate has significantly reduced 
the large theoretical uncertainty present in the previous LO calculation. From the quark- 
level b — > S7 decay rate it is possible to infer the B meson inclusive branching ratio BR 7 = 



BR(i? — > Xg'y) by including the small (percent-level) non-perturbative l/m& ]I3 and l/m c 
fl9|j , |p0|| corrections. The result, the most recent theoretical prediction for BR^ LO in the SM, 
is 
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10 4 BR^ LU = 3.48 ± 0.31 (31) 

where the theoretical uncertainty, ~ 10%, is less than half the previous LO uncertainty. It 
includes the uncertainty in the SM input parameters m t , a s (Mz), a em , m c /mb, rrib, CKM 
angles and the residual scale dependence uncertainties. Nevertheless the error is dominated 



by the uncertainty on the SM input parameters, as claimed in Ref. p2 |. Comparing with 
the experimental result in Eq. (|2"9"|), the SM prediction is higher than the observed branching 
ratio, although the disagreement is less than two experimental standard deviations. 

To compute the deviations expected in the MGM, we look closer at the theoretical calcu- 
lation. The non-leptonic low energy effective Hamiltonian relevant for the b — > S7 decay 
is 

rNL . 4GV 



^eff 



-^WEWftOi*) (32) 
V2 i= i 



where are the CKM matrix elements, Cj(/^) the Wilson coefficients and Qi(fib) the cor- 
responding operators evaluated at the renormalization scale fib — 0(rrib). Only the magnetic 
and chromomagnetic dipole operators, Qr and Q$ respectively, are significantly affected by 
new MGM amplitudes; their expressions are given by 

= T7HI m b(s L a^b R )F^ (33) 
lo7r z 

Qs = -^- 2 m b {s L a^T a b R )Gl v (34) 

where T a are the SU(3) co ior generators and and G" are respectively the electromagnetic 
and SU(3) co i or field strengths with e and g s the corresponding coupling constants. The remain- 
ing operators in the complete basis can be found in Ref. |H] or [f^| . The Wilson coefficients 
Ci(n), which satisfy the RG equation 

/^Ci(^) = CjWijiin), (35) 
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can be expanded, together with the anomalous dimension matrix 7^, in powers of a s as follows: 
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47T 
t(0). 



Then the solution in the MS renormalization scheme for = is given by 



(36) 
(37) 



Cffo) = ^Cr(M w ) + ^(^-r / i)CT(M w ,) + ^M ai 

^ i=i 

Z297664 is 7164416 u 256868 37 6698884 39 \ , m . 
+ I -r-^T^rV 23 ocr ^ n ^ cr ?7 23 + ; /looo ?7 23 t^t^tV 23 Q (Mv) 



(38) 



V 14283 



357075 



14283 

8 



357075 



+ 



• J> ' 208 - ^M) cf (M w ) + J] (e^fo) + /, + 



4761 



(39) 



where 77 = a s (Mw) / <x s (fJ>b) > x 4 = m 2 /M^, and a s is calculated at NLO. The SM contribution 
to C?s(Mw) together with the real numbers fi, hi, and the function E(x) can be 
found in Ref. El. 



Finally the branching ratio BR 7 , conventionally normalized to the semileptonic branching 
ratio BR(5 -> X c eis) = (10.4 ± 0.4)% |23|, is given by |2T 



BR 



NLO 



BR(5 -> X c ei/) 



|V^14| 2 6a em 
|V; b | 2 ixg(z)k(z) 



3 7T 



47T 



V i=i 



/ \ 1 (0) , m fc 



A = (e- Qs ^ los5 ( 7+21os5 )/ 37r -l)|Cj 0) (/i fe )p 



+ 



71 



(40) 



where z = m^/m 2 and in the last sum i < j. The expressions for all LO Wilson coefficients 
Cj; \nb) together with the functions g(z), k(z), r^z) and fij(S) can be found in Ref. []2T|. The 
term 5 np , of order a few percent, includes the non-perturbative l/m& 18fl , and l/m c fL9 |, pO 
corrections. 



Some comments about this result are in order: 



1. With NLO accuracy the inclusive b — > 57 decay rate is given by the sum of the total 
rate for b — > 57 decay and its gluon Bremsstrahlung corrections b — > s^yg, where in the 
former the 0(a s ) radiative corrections to the matrix element are included and in the 
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latter an explicit lower cut on the photon energy resolution J5 7 > (1 — 5) 777.5/2 is made. 
In the sum of these two contributions the only remaining infrared cutoff is the logarithm 
of 5 in the A term of Eq. (gUp. 

2. The LO result can easily be recovered by taking the limit of a s — > inside the expression 
for the D and A terms. 

3. To retain a strictly NLO result, the terms proportional to a 2 in the \D\ 2 expression 
should be discarded. 



Since the new physics scale is above the electroweak scale, the new contributions to the 
b — > 57 decay will affect only the SM Wilson coefficients at the M\y scale. Therefore the non- 
perturbative resummation of the large logarithms from My/ to fib can be taken completely 
from the existing SM calculation. Of course this would not hold if there were new operators, 
generated at the electroweak scale, which would significantly correct C?([Ab) via mixing with 
Qj. But in the MGM this does not occur: the dominant contributions to the h —> 37 decay 
enter only in the Wilson coefficients C 7 (M W ) and C 8 (M W ). We define fractional deviations 
i? 7 8 from the SM amplitudes: 

cfl{M w ) ee C\f U {M w ) (1 + R r , s ) (41) 

where Cf^ M (My/) represent the LO SM contribution. Inserting these definitions into the 
BR 7 formula in Eq. ( |4*0|) yields a general parametrization of the branching ratio in terms of 
the new contributions UTTl : 



10 4 BR^ LO = (3.48 ± 0.31) (l + 0.622f2 7 + 0.090^ 



+ 0.066i? 8 + 0.019R 7 R S + 0.002i?y (42) 

where the following central values are used: mf ole ~ m^ s (m^) ~ 174 GeV, m^ ole = 4.8 GeV, 
m poie _ 13 QeV, fib = rrib, a s {mz) = 0.118, a~^{mz) = 128, sin 2 ^ = 0.23 and a photon 
energy resolution corresponding to 5 = 0.9 is assumed. 

What are the contributions to Rr^ m the MGM model? In general in the MSSM the 
diagrams which contribute to the b — > S7 amplitude, in addition to the SM, are given by a 
chargino, neutralino, gluino and charged Higgs exchange. 



The charged Higgs-exchange amplitude (*4h-) is well-known to be large p^| , |L5[], [{17 
More recently it was pointed out that for large tan (3 also the chargino and gluino exchange 
can give sizeable contributions to the b — > 57 amplitude fl5 |. In particular, the leading large 



tan/3 contributions to the amplitude arise from exchanges of a pure higgsino (A.^-), a mixed 
wino-higgsino (A^-), or a gluino (Ag), in descending order of importance. Thus 

Asm 
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where the amplitude *4sm correspond to an exchange of a W boson. We stress here that the 
same mechanism of tan/3 enhancement is present also in the 5m b /m b radiative corrections to 
the bottom quark masses [15]. Indeed the diagrams contributing to the b —■ sj amplitude are 
the same as for 5m b /m b , but where a photon line is attached in all possible ways. Also, as 
explained above, the MSSM contributions are proportional to the relative sign of /x, My 2 and 
tan fi. In the MGM this sign is calculable, and therefore a unique prediction of b —■ sj decay 
rate in terms of only the messenger and gaugino mass can be given flO| . JlTJ . 



In the MSSM the complete and exact analytical results for the b —■ sj and b — > sg total 
amplitudes can be found in Ref. [H], in terms of the various mass eigenvalues. However, 
since in the MGM the squarks are quite heavy, a mass insertion approximation will suffice 
TO]. Here we diagonalize the chargino mass matrix exactly, since we are interested in light 



charginos. The result for the partial amplitudes is, up to an overall normalization [\ I 



A 



SRI 



A u - 



A 



Wh- 



V ts XtW 



-VtsHXtH 



'-Fx{xtw) + F 2 (x tw ) 
;F 3 (x tH ) + F A (x tH ) 



-V ts r b tanf3 



t^tR 



1 1 
tL tR 



Mr 



X 



hit) 



Ml 

hi 



Mi 



r b tan fi- 



lm 



m~ m~. T 



M r - M 'r 



h 



a. 



-r b t&nfi 
9 a w 



rrvLmj ~ 

w t L c L 



fj,M 3 



rrr- 



F z \(xm- 



;iqj 



F'{ 



X 



im? mz 

I tL t R \ 



h\ii X h 2 t) 



M? 

h. 



Ml 



(44) 
(45) 

(46) 
(47) 
(48) 



where r b = 1/(1 + 5m b /m b ) accounts for the mass corrections to m b fllQ |. The functions Fi 
can be found in Ref. HT4|, while the new functions are defined as follows: 



F(x 1 ,x 2 ) 
/(*) 



f(xx) - f(x 2 ), F'(xi,x 2 ) = Xi/(a;i) - x 2 f(x 2 ), 
d , „ 2 „ s „ _ Id 2 



-(xF 3 + ^xF 4 ), F gl (x) = §^(^ 2 ^). 



(49) 
(50) 



We also use: 



XtW 



X hiq L 



m. 



m 



2 ' 

w 
Ml 



XtH 



m 



H- 



X hit 



m L m c L \ 



XM 3 q 



2 



Ml 

h. 



\ m tr m t f 



(51) 



(52) 



where Mz- is the zth chargino mass eigenvalue and m„ is the average squark mass. The left- 

i 

right stop mass insertion is $ — Tn t A t (> 0), while the left-left stop-scharm mass insertion 
is given by ]T0J 

A Mn(^). (53) 



VTLr - ~ 
tL^L 



+V ts mj 



4tt 2 



mz 
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The contribution to Rg can be obtained from Eq. (13) after the following replacement in 
the amplitudes of Eqs. (Pj-pD : 



\f x +F 2 ^F x , ?-F 3 + F 4 ^F 3 , /^A(xF 4 ) (54) 
6 6 ax 

Fgl{x) ~" Yq^ {3x2f3 + \ x2Fa) - (55) 

At present in the MSSM the super symmetric corrections to C-i(Mw) and C%{Mw) are 
known only at LO; we will return to comment on the resultant uncertainties when we present 
the numerical predictions for BR 7 decay, in terms of the SUSY breaking scale and the mes- 
senger mass scale M M , in Sec. 6. 



4 b -> s£ + £~ 

In this section we consider the semileptonic flavor-changing decays b — ► se + e~, b — > sfj, + ^i~ 
and b — ► sr + r~. These processes are interesting for several reasons: first, they are sensitive 
to new physics with large tan/? for the same reason that b — > 37 is sensitive, and in fact 
because of the same operator Q 7 (see also Ref. p0| ); second, they involve other operators as 
well, and so can serve as complementary tests of the model [25|; and third, since there are 



several measurable quantities (branching ratios and asymmetries in several systems), they can 
potentially yield much more data. At present these decays are known at the NLO logarithmic 



accuracy for the SM p7 |, [28|. The l/m& non-perturbative contributions are small and well 
under control, except near the end-point of the dilepton mass spectrum pEf . Recently it 
was pointed out that other non-perturbative power corrections 0(h^ CD /m 2 c ) could affect the 



amplitude, but in practice they are very small [pO] . Moreover there are other non-perturbative 
effects due to the resonance regions in the dilepton invariant mass distribution. But these do 
not have to be taken into account if certain regions are excluded from the integration region. 

Unfortunately, the expected branching ratio for these decays is too small to be observed in 
present experiments. The SM predictions for the non-resonant decay rates BR(6 — > s£ + £~~) = 
BR« are given by p9| , |IT| 



BR ee 


~ 7.0 x i(r b 




~ 4.5 x 


BR rT 


~ 2.0 x 10" 7 



(56) 
(57) 
(58) 

at present, experiments only place some upper bounds on BR ee and BR W , 0(few x 10 -5 ), or 
about an order of magnitude above the SM |32| . However an increase of two to four orders of 
magnitude in the experimental sensitivity is expected in the near future at the Tevatron |54" 



With such statistics it will be possible to measure the branching ratio and the energy or 
angular distributions, even for the more rare r decay mode. 
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We now turn to our prediction for the branching ratios and asymmetries. The effective 
Hamiltonian H e Q relevant for the b — > s£ + £~ decay is: 

Gi 
V2' 



H, 



cff 



H^ L - ^v*v tb {CMQM + CMQM} 



(59) 



where is given in Eq. ( j3^) with ~ 0(m&), the magnetic-dipole operator Qi is defined 
in Eq. (^) and the semileptonic operators Qg and Qio are given by 

(60) 
(61) 



Qg = (s^(l 
Qio = (s^(l 
The Wilson coefficients Cg(fi b ) and Ci (fi b ) at NLO are: 



7b)6) (/yi 
7s)6) (h^i 



a 



C g (ji b ) = —(P Q (ji b ,ri)+P E E) 

Z7T 



Cg(M] 



w 



where the initial conditions Cg j io(Mw 

Cg(M w ) = 



in the HV scheme are given by 
1 



a 
2^ 



sin 2 fl 



-Y-4Z 



w 



C 10 (M } 



a 



w 



-Y. 



(62) 

(63) 
(64) 



2tc sin 9w 

The terms Y, Z and E contain the results of the one-loop integration for the 7, Z penguin 
and box diagrams at the electroweak scale. The SM expressions for Y, Z, E, respectively the 
functions Y SM (x t ), Z SM (x t ) and E SM (x t ) (x t = m 2 t /M^), can be found in Refs. [Z7|, The 
function Pq includes the resummation of the large logarithms at the NLO and it is scheme 
dependent. Finally, in the limit 77 — > 1, Pe —>■ and Po — > 0, so that the initial condition 
Cg(fi) — > Cg(M w ) is recovered. The expressions for P and Pe, in both HV and NDR scheme 
for the former, can be found in Ref. [^8[ . 



Note that Cio(//&) does not depend on the renormalization scale \i b [see Eq. flBTDl. This 
is a consequence of the P and C symmetry of the electromagnetic current which forbids any 
mixing of Qio with Qg or any other operator of the non-leptonic effective Hamiltonian H^. 

At NLO accuracy the total amplitude for b — > s£ + £~ decay can be expressed as a function 



a. 



of Cf(fi b ) = C 9 (fi b ) + C 9 MB (/i fe ), CH/ife) and C 10 (M W ) where C 9 ME (n b ) contains the 0{ 
one- loop matrix elements needed to cancel the leading \i b dependence in C9 (//&). Moreover 
Cg S is scheme independent since C g ME cancels exactly the scheme dependence of Cg(/x&) p8fl . 
Note that at NLO accuracy, C>?(pb) enter only at the leading order. 



Finally the differential decay rate for b 
d 2 T Pf 



s£ + £ is given by 



dy+dy. 



(~<2 

F 5 I t /* T / |2 2 



+ 
+ 



256tt 5 



°9 



Oft 



K 



7,7 



Cf(/i 6 )Re 



i^ 9,9 + C 2 (M W )K 10,10 

0&) 



C 10 (M W ) (C^\fi b )K 7 ' 10 + Re [Cfifj, 



Ceff 
9 



K 7 ' 9 

K 9,10 



(65) 
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Note that (//&) is real [see Eq. (|38|)1 as is Cio(Mw), while Cg S (nb) is complex — its imag- 
inary part originates in the one-loop matrix elements in Cg IE (fib)- The kinematic variables 
are y± = 2E±/mt, where E± are the lepton energies measured in the b rest frame. We prefer 
to use the following combinations: 



y+ + V- - 1 



(66) 



y 



cos 9 



y+ - y- 

where y max — (1 — s)J 1 — Ax^/s, xi = m\jm\ and 6 is the angle between the b and the £ + in 
the £ + £~ rest frame. Including only the lepton mass corrections, we find 



#9,9 


1 

2 


1 - s 2 - y 2 + 4x e (l - s) 


(67) 


#-10,10 


1 
2 


\-s 2 -y 2 - 4xt(l - sj 


(68) 


#7,7 


2 

s 


'l-s 2 + f + ^{ 

s 


1-5) 


(69) 


#7,9 


= 4 


l — s+ —(1 - s ) 




(70) 


#7,10 


= # 






(71) 


#9,10 


= 2sy. 




(72) 



where the massless case is 



The results of Eqs. (0), (|67|-[72]) are in agreement with Ref. 
considered. Moreover, after integrating Eq. flBSp with respect to y, our result agrees with that 



of Ref. [31], including the mass corrections. 



The 1/s terms present in Eqs. ( ]69|J70| ) are due to the photon propagator which connects 
the quark matrix element of the magnetic dipole operator Q-j to the electromagnetic lepton 
current. As observed in Ref. [2~9"| the last term of Eq. ( p5|) proportional to the lepton mass 
(xi) can not be completely neglected even in the light lepton case. In particular it gives a 
finite contribution to the total branching ratio (roughly a few per cent) even in the m 2 — > 
limit. Indeed the 1/s 2 term in Eq. (|69]) generates a pole 0(l/mf) in the integrated branching 
ratio which is cancelled by the 0(mf) terms present in the numerator. 

The total branching ratio for a particular I is obtained integrating d 2 T 'a/ 'dy + dy_ over the 
entire kinematically-allowed range for that I. The forward-backward asymmetry A« (which 
is the same as the energy asymmetry) is proportional to the integral over the range y_ > y + 
(i.e., cos# > 0) minus the integral over the range y_ < y + (cos 6 < 0). The allowed kinematic 
limits are -y max < y < +y max and Ax e < s < 1. 

To avoid intermediate charmonium resonances and non-perturbative phenomena near the 
end point, we exclude the same ranges of s specified in Ref. J29|. In particular, the integrating 
regions of s in the I = e, \x case are 



(m\s) E {Am 2 , (2.9 GeV) 2 } U {(3.3 GeV) 2 , (3.6 GeV) 2 } U 
{(3.8 GeV) 2 , (4.6 GeV) 2 } 



(73) 
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Figure 1: Dalitz plots for the integrating regions of b — > s£ + £ decays corresponding to e, \x and 
r decay modes from the left to the right respectively. 



and for I 



(mis) e [Am 2 T , (3.6 GeV) 2 } U {(3.8 GeV) 2 , (4.6 GeV) 2 } 



(74) 



The Dalitz plots for the decays b — > s£ + £~ are shown in Fig. |I|. The dimensionless lepton 
energies y + and y~ run along the horizontal and vertical axes, respectively. The rotated axes 
correspond to y (pointing right and down) and s (pointing right and up). Only the regions 
of s included in the integration [see Eqs. ( |?3| , |74"D 1 have been shaded. The branching ratios are 
obtained by adding the integrals over the (colored) dark- and light-shaded areas, while the 
asymmetries require subtracting the integral over the light-shaded region from the integral 
over the dark-shaded region. 



As usual, we normalize the branching ratio to the semileptonic decay rate b 
order to eliminate the large dependence on the b quark mass: T« — > Yu = ^u/^b- 



— > ceu e m 
. c , hence 



(BR, 



b^c 



ds2 



aim 



\ dsdy ) 



(75) 



symm 



N(y_ > y + ) - N(y+ > y_) 
iV(y_ > y+) + N(y + > y_) 



ds2 



Vmax 



dy 



f d 2 f u \ 
\ dsdy J 



(76) 



antisymm 



in which the symmetric (antisymmetric) rate must include only the kinematic coefficients 
symmetric (antisymmetric) in the exchange y — ► —y. 

We have also calculated some polarization asymmetries, namely the rate when one of the 
lepton's spin is parallel to its momentum versus the rate when that spin is antiparallel to the 
momentum. We have considered both the integral of the asymmetry in the differential rates, 
and the asymmetry in the integrated rates, and also looked at the three possible polarization 



directions, using the work of ref. |33| . 
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Let now discuss the impact of the MGM model in the b — *> s£ + £ decays. 

The contribution of the MSSM to C$ } \a(Mw) is given by the 7 and Z penguin diagrams 
proportional to the form factor 7^, and by the box diagrams. Their contribution are incor- 
porated into the Y, Z and E expressions. Inside the penguin and box diagrams can run 
the charged Higgs, charginos, gluinos and neutralinos. Since Pe is two orders of magnitude 
smaller than P , and the SUSY contribution to E should be smaller or comparable to the SM 
one, we will neglect the E contribution in our analysis. The complete analytical results for 
Y SUSY and Z SUSY can be found in Refs. 0, |2|. 



In the MGM model we have computed the SUSY contribution due to chargino, gluino and 
neutralino exchange in the 7 and Z penguin diagrams and box diagrams and we find that 
Y SUSY and Z SUSY are very small comparing to the leading SM contributions, roughly a few 
percent. (Note that the squarks are quite heavy, but there is no compensating large tan/3 
enhancement.) Also the charged-Higgs contribution to Y and Z is small, being suppressed by 
at least 0{m s /m b ) and in practice it is of order of a few percent in both the penguin and box 
diagrams. 

New chirality- violating four-fermion operators could be generated at the squark mass with 
their chiral suppression partially compensated for by tan/3. However we found that their 
contribution is small compared to the leading SM contribution to C 9tW (M w ), at least for 
large squark masses 0(0.5-1 TeV). 

Thus we find that the only relevant effect of the MGM on the b — > s£ + £~ decays can be 
parametrized by R-j in Eq. (|43|). In particular for the BR« and An we find [11] 



10 7 BR ee = 68.5 + 22.4^7 + 6.1^ (77) 

A ee = (4.52-3.01i? 7 )/(10 7 BR ee ) (78) 

10 7 BR W = 44.64 + 2.46^7 + 1.81^ (79) 

A wt = (4.68- 2.90i2 7 )/(10 7 BR W4 ) (80) 

10 7 BR TT = 2.013 - 0.20LR 7 + 0.009^ (81) 

A rr = (0.434 - 0.042 J R 7 )/(10 7 BR TT ). (82) 



The coefficients of these Rj polynomials were obtained using the same central values for the 
SM parameters used in the b — ► 57 decay. We have checked the sensitivity of these predictions 
to the various input parameters and to the scale fif,. The main sources of uncertainty are 
m t , a s {m z ) and the residual ^ dependence. We find that these uncertainties affect the 
normalizations of the BR« by up to (9(10%). However the uncertainties on the functional 
dependence on R7 are much less, at the percent level. There is a small sensitivity to Rs (which 
nonetheless can be included) . The normalizations of A ee and are even more sensitive than 
the branching ratios to fib, but an improved SM NLO calculation can reduce the uncertainties 
to less than ~ 10%. On the contrary the r forward-backward asymmetry A rr has a weaker 
functional dependence on Rj then the other decay modes. Moreover we find that the various 
integrated polarization asymmetries are not particularly sensitive - at most at the few percent 
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level - to i?7, so they would not serve to test this model and we do not include them in our 
final results. 

Finally we would like to stress that, when the present residual \ih dependence is reduced 
by improving the SM NLO calculation, and when more precise measurements of rrit and 
oi s {mz) are available, then the MGM model predictions for BR« and A« can be significantly 
sharpened. 



5 g,~2 

The anomalous magnetic moment of the muon a M = (g M — 2)/2 is one of the most important 
high precision experiments, providing extremely precise tests of QED and electroweak interac- 
tions as well as strong contraints on new physics models. The theoretical prediction is known 
to 0(a 5 ) in QED and recently also the two-loop electroweak radiative corrections have 
been included P6| |. The agreement between the theoretical prediction and the experimental 
result is impressive - but not exact. The current average experimental value of a^ xp is |3T 



10 10 < xp 



11 659 000) = 230 ±84 



13) 



The E821 experiment at Brookhaven National Laboratory (BNL) is expected to measure a M 
to within ±4 x 10~ 10 , or perhaps even ±1 — 2 x 10~ 10 []3~8|| . The most up to date theoretical 



prediction in the SM is given by the sum of the following contributions [39 ] : 



10 



1CKSM 



OED + fl EW + a H(l) (vac pol) + m 2)(vac pd) + a H, x 7) 



*4) 



where 



-.QED 



,EW 



7^(2), 



vac pol) 
vac pol) 

<(77) 



11 658 470.57 (0.19) 
15.1 (0.4) 
701.1 (9.4) 
-10.1 (0.6) 
-7.92 (1.54). 



(85) 
(86) 
(87) 
(88) 
(89) 



the most recent pure QED contributions at order a 5 ; aj^ w includes the two- 
loop |36| electroweak corrections; the a^ 1,2 ^(vac pol) include respectively the one-loop |4T 



qed i nc i uc i es 



and two-loop J|2 hadron vacuum polarization corrections extracted with dispersion relations 
from the measurement of the e + e _ — > hadrons cross section; and (77) contains the hadronic 
"light-by-light scattering" contribution f43|. 



The hadronic contribution to the vacuum polarization is the largest source of error, re- 
flecting the large experimental uncertainty of the data. However with the new experiments of 
BEPC at Beijing, DA$NE at Frascati and VEPP-2M at Novosibirsk the theoretical error in 

(vac pol) will be significantly reduced [ |3| . The hadronic light-by-light contribution arises 
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from a different class of diagrams which cannot be directly extracted from experimental data, 
and therefore a somewhat model-dependent calculation is used. 

Finally, by combining the results of Eqs. (|34|-|S9l) and summing the errors in quadrature 
we arrive at the following SM prediction: 

10 10 aJ M - (11 659 000) = 170 ± 10 . (90) 

The improved experimental precision in the next experiment at BNL and the reduction of the 
theoretical uncertainties in the hadronic vacuum polarization will allow a direct test of the 
electroweak corrections. 

The expected deviations from a^ M induced by the MSSM, and the degree to which the 
MSSM can be constrained by improved experimental measurements, has been studied since 
the early days of supersymmetric phenomenology |[45||. More recently the impact of large 



tan/? was considered, first in the context of gravity- mediated SUSY breaking [46], next in 



the general MSSM [J7|, and finally specifically in gauge- mediated scenarios |48|| . However, 
in this last work, the crucial parameters /x, tan/3 and the sign of the new amplitudes were 
not predicted by the model. In contrast, in the MGM model we consider all three of these 
crucial parameters are predicted along with the rest of the MSSM spectrum in terms of our 
two fundamental parameters, allowing us to correlate our prediction with the predictions 
of b — > S7 and b — > s£ + £~. 

The MSSM contributes to mainly via magnetic-dipole penguin diagrams, with an ex- 
change of a chargino or a neutralino in the loop. We have calculated these contributions and 



our analytical results are in agreement, in magnitude and sign, with the results of Ref. |47 
(after accounting for the erratum). 

In the MGM model the chargino amplitudes dominate the neutralino ones, and within 
the chargino contribution the terms proportional to the chargino mass completely dominate. 
Using also <C vjj yields |pT |, for the deviation from the standard model, 



where 



A< OM a £ K (91) 

" An mg(A/?_-M?_) v 1 



F ll (x 1 ,x 2 ) = Ufa) - fpfa), (92) 
3 - Ax + x 2 + 2 In x 



/«(*) 



3(1 -x) 



and x~ hi ~ = M^-jm-. For our numerical results we will use the complete expressions for 



MGM 



6 MGM model predictions 

As explained in Sec. 2, when the self-consistent equation Eq. (|2~4] ) for tan/3 is numerically 
solved, all of the relevant SUSY spectrum and physical quantities can be predicted as a func- 
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tion of only two fundamental parameters: the effective SUSY-breaking scale A or equivalently 
the weak gaugino mass M 2 , and the logarithm of the common messenger mass Mm- In prin- 
ciple these results will depend on the doublet and triplet messenger mass splitting, but the 
dependence is very weak for the expected small splitting. In all our computations we will 
assume that Mm 3 = 1.3Mjw 2 . 

We consider the range 2A < Mm < 10 4 A. Messenger masses much nearer to the absolute 
lower bound A (below which their scalar components would develop VEVs) are not of much 



interest: they are either already ruled out by vacuum stability constraints |L0| or require 
extremely heavy superpartners and thus mimic the standard model. Messenger masses much 
above (9(10 4 A) may also destabilize the vacuum or run into cosmological difficulties with a 



heavy gravitino (110, BUf and references therein), though this is by no means an airtight bound. 



Within this range of M M , the one-step solution we have used for the RG equations is 
an adequate approximation. In the low M M region, tan/5 is ~ 50 and there are strong 
cancellations in the B term and in the b — > S7 amplitudes, so an accurate solution is required. 



But as shown in Ref . |1(J , the one-step approximation in this case is as accurate as keeping only 
the leading threshold corrections, and both of these approximations are entirely sufficient for 
our purposes. For the highest values of Mm we consider, the one-step approximation is not as 
accurate, but the cancellations in B and in b — ► 57 are not as severe, so the reduced accuracy 
is once again sufficient. In the end, we estimate that the uncertainties in our predictions of 
BR 7 and comparable to the corresponding SM ones. 

In Figs. |2| and ||| we present our MGM model predictions as a function of Mm and the weak 
gaugino mass M 2 for = 1 or 2 messenger families, respectively. Different curves are given 
for different messenger masses M M , and heavy (red) dots indicate discrete values of M 2 . The 
SM-predicted and the measured central values are shown as circled dots. The dark and light 
green areas surrounding the experimental central value indicate the regions allowed at 1 a and 
95% CL, respectively. The dark ellipse surrounding the central value of the SM predictions 
indicates the expected theoretical uncertainty, which should also be applied to each point on 
the MGM prediction curves. Finally, the small dashed ellipse centered (for the present) on 
the experimental central value indicates the 1 a uncertainties expected in the near future from 
BNL and CLEO. 

The heavy (continuous and dashed) red curves are our predictions for a M along the vertical 
axis, and either BR 7 or any of the various B — > X s l + l~ branching ratios and asymmetries. 
All the latter are correlated in their dependence on the single parameter R7 which measures 
the deviation from the SM of the amplitude of the Q7 operator. At present BR 7 is by far 
the best measure of R7, which is why we emphasize it in our graphs. (The various quantities 
along the horizontal axes depend slightly and differently on the R$ parameter, so to present 
them in the same graph we have chosen a specific value of Rg = 0. But we estimate that this 
approximation affects the R? predictions by an uncertainty that is well within the SM one 
shown in the figure.) The curves become dashed for those ranges of Mm and M 2 which are 
in fact ruled out by vacuum stability constraints ||10|| , primarily due to the stau developing a 
charge-breaking VEV. 
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Figure 2: The N = 1 MGM predictions (as heavy red curves) of BR 7 (xlO 4 ) and lO 10 ^ - 2) - 
1165900 as functions of the wino mass M2, for various messenger masses Mm and with heavy red 
dots at discrete M2 values. The experimental and SM central values, indicated with circled dots, are 
surrounded by the corresponding uncertainty ellipses. The extra horizontal axes are the correlated 
predictions for R7, BR« (xlO 7 ) and A«. 

Several features of the MGM model predictions stand out from these figures: 

1. When the messengers are relatively heavy, the MGM predictions are in better agreement 
with experiment than those of the SM (except for extremely light superpartners): when 
the superpartners are lowered (to their natural values, around the electroweak scale) and 
the MGM departs from the SM, a M increases and BR 7 decreases towards the experimental 
central value. In particular, the sign of this effect is a prediction of the MGM. 
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Figure 3: The MGM predictions as in Fig. | but for N = 2. 

2. When the messengers are quite light, a fortuitous cancellation between the chargino and 
charged-Higgs amplitudes in R7 prevents BR 7 and the various B — > X s l + l~ branching 
ratios and asymmetries from deviating significantly from the SM, although a M is still 
raised as before. 

3. The expected reduction in experimental uncertainty will tremendously sharpen the com- 
parison with the SM and MGM predictions, possibly deciding between the two and 
definitely strongly constraining the latter. 

4. At present, the (absolute) bounds vaccum instability bounds on Mm and M 2 are stronger 
then the (95% CL) experimental bounds due to BR^, and a p . 
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As a function of N the sfermion masses scale like 1/yN relative to the gauginos [see 
Eqs. (0) - ©]- Raising N for fixed gaugino masses thus lowers the sfermions and amplifies the 
MGM effects. While this is the general trend, in practice there are some complications because 
of cancellations and the variation in tan f3. By comparing Figs. |^, |] throughout the allowed 
region for fixed values of Mm and M 2 , we see that when N is increased, R-j (or equivalently 
BR 7 ) is decreased and a M is increased. Raising iV is roughly equivalent to lowering M 2 while 
raising Mm- 

From the horizontal axes below those of BR 7 and R-j in Figs. 0, ^, our predictions for BR ee , 
BR TT , A ee and A MM may be read off. We do not give the results for BR W and A TT because 
they are not very sensitive to R-j. Evidently the MGM model can produce large deviations 
from the SM in A ee and A^, and somewhat smaller ones in BR ee and BR^. Note however 
that the b — > s£ + £~ quantities have their own theoretical uncertainties beyond those indicated 
in the figures, mainly from the residual dependence, as explained in Sec. 4. These could 
be reduced by improving the SM NLO (in a s ) calculations, at which point the MGM model 
would predict the branching ratios and asymmetries more precisely. We expect that at least 
some of these predictions could be tested at hadronic colliders in the next few years. 

7 Conclusions 

Within the minimal gauge-mediated SUSY-breaking model, which can naturally generate a 
large tan j3, we have computed the inclusive branching ratio B — ► A s 7, the inclusive branching 
ratios and asymmetries for B — ► X s £ + £~ (with £ = e, r), and the anomalous magnetic 
moment of the muon = (g^ — 2)/2. In particular we included the complete next-to-leading 
order accuracy in the strong coupling for the B decays. 

The MGM model, described in Sec. 2, has a naturally high tan/? signature provided 
that the tree-level soft-breaking B parameter, which couples the Higgs doublets in the scalar 
potential, vanishes. Moreover this model is highly predictive: all the SUSY soft-breaking 
terms, as well as tan j3 and the physical relative sign (often called the sign of the \i parameter), 
can be predicted in terms of only two fundamental parameters: the SUSY breaking scale (or 
equivalently the wino mass), and the logarithm of a common messenger mass. Therefore the 
MGM predictions for these processes are strongly correlated. 

The correlated predictions are shown in Figs. ^| and |3|. Note that the correlation between 
a M and BR 7 depends on the model parameters, Mm or M 2 , whereas the various b-quark-related 
quantities are all completely correlated through their dependence on the single quantity R 7 
independent of M 2 and M M (a quite general feature of the MSSM, as discussed in Ref. Q). 
We have gone slightly beyond the minimal model by including the results for N = 2 messenger 
families. We included the results of the numerical predictions for N = 1 and N = 2 messenger 
families in Figs. |2] and |3] respectively. The general trend of the MGM predictions emerging 
from these results can be summarized as follows: 
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• When the messenger scale is a few orders of magnitude above A ~ (9(100 TeV), and 
the superpartners are near their natural electroweak scale, the MGM predicts a higher 
anomalous magnetic moment and lower b — > 57 branching ratio than the SM, in agree- 
ment with present measurements. 

• When the messengers are not far above A, the superpartners must be significantly higher 
than the electroweak scale, and (due to a cancellation) the b — > 57 rate is largely 
unperturbed from its SM value, though the anomalous magnetic moment is raised. 

• Increasing the number of messenger families amplifies these predictions somewhat: it is 
similar to lowering the superpartner scale and raising the messenger scale. 

• At present, the allowed range of M 2 and Mm values is more constrained by theoretical 
vacuum- stability constraints than by experimental data on BR 7 and a^, though that 
data does mildly favor the MGM model over the standard one. However, the BNL 
experiment now in progress, and the much-awaited CLEO analysis of their BR 7 mea- 
surements, will dramatically constrain the MGM and allow a much more convincing 
discrimination between the two models. Measurements of several B — > X s £ + £~ branch- 
ing ratios and asymmetries at hadron colliders should further sharpen the b-quark side 
of this picture in the next few years. 



Note added: After this work was completed, new and significant theoretical results on 
radiative b decay have appeared, one [5(| taking into account the proper electroweak radiative 
corrections and the other [|51|] convincingly criticizing the current method of extracting the 
inclusive rate for b — > 57 from the currently published CLEO data. The former slightly lowers 
the SM prediction, while the latter argues that theoretical uncertainties have so far been 
underestimated, and that only with more new data on the spectrum or with an improved mi, 
value from upsilon spectroscopy can the uncertainties be reduced. Thus, though at present it 
is difficult to sharply compare the predicted inclusive decay rate with experiment, we expect 
that in the very near future such a comparison will be possible. 
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